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We present a new calculation of the D —> -k,Iv semileptonic form factor f+^^iq 2 ) at q 2 = 
based on HISQ charm and light valence quarks on MILC Nf =2 + 1 lattices. Using methods 
developed recently for HPQCD's study of D K, lv decays, we find /+~* 7r (0) = 0.666(29). This 
signifies a better than factor of two improvement in errors for this quantity compared to previous 
calculations. Combining the new result with CLEO-c branching fraction data, we extract the CKM 
matrix element \V c d\ = 0.225(6) e:E p. (lQ)i a t. , where the first error comes from experiment and the 
second from theory. With a total error of ~ 5.3% the accuracy of direct determination of \V c d\ from 
D semileptonic decays has become comparable to (and in good agreement with) that from neutrino 
scattering. We also check for second row unitarity using this new \V c d\, HPQCD's earlier |V C3 | and 
\Vcb\ from the Fermilab Lattice & MILC collaborations. We find \V cd \ 2 + \V CS \ 2 + \V cb \ 2 = 0.976(50), 
improving on the current PDG2010 value. 

PACS numbers: 12.38.Gc, 13.20.Fc 



I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix pro- 
vides particle physicists with a wealth of opportunities 
to carry out precision tests of the Standard Model (SM) 
and look for New Physics. On the one hand each ma- 
trix element can be determined in several ways, em- 
ploying different experimental and theory inputs, and 
the results compared with each other. Three gener- 
ation unitarity can also be examined to see how well 
Vckm x V(j KM = I is satisfied. This leads to tests 
such as 1 st , 2 nd or 3 rd row/column unitarity. It also 
gives rise to the important "Unitarity Triangle" relation 
V ud * V* b + V cd * V* b + V td * V* b = 0. Consistency checks of 
the sides and angles of the Unitarity Triangle (UT) have 
been the focus of much of the experimental and theo- 
retical effort in Flavor Physics in recent years. Lattice 
QCD is playing an increasingly important role in CKM 
physics [1]. For instance, lattice calculations of the Kaon 
semileptonic form factor /^"^(O) [2, y] and the decay 
constant fx (or Jk/ fir) 043 have contributed to preci- 
sion determinations of \V US \ and 1 st row unitarity tests 
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The HPQCD collaboration recently published a new 
lattice calculation of the D — > KAv semileptonic decay 
form factor f+^ K (q 2 ) at q 2 = [10| which significantly 
reduced the error on this quantity compared to previous 
theory results. This led to a very precise determination of 
the CKM matrix element |V^ S |. Features in the HPQCD 
work that made this improvement possible include the 
use of a relativistic quark action, the "Highly Improved 
Staggered Quark" (HISQ) action [ll[, to simulate both 
light and charm quarks, better data analysis tools and 
a new method for carrying out chiral/continuum extrap- 
olations. Capitalizing on these developments, we turn 
here in this article to D — > 7r, lv semileptonic decays. We 



focus on extracting the CKM matrix element \V c d\ by 
combining theory results for fF^ n (0) with experimental 
input from CLEO-c [Hj]. The first unquenched lattice 
studies of D semileptonic decays were carried out several 
years ago by the Fermilab Lattice & MILC collaborations 
In that work lattice gauge theory was able to pre- 
dict the shape of f+{q 2 ) prior to subsequent confirmation 
by experiment. The theory errors for fS^ w (0) in (l3j 
were ~10%, and this has remained the dominant error 
in determinations of \V c d\ from D semileptonic decays. 
More accurate determinations have come from neutrino 
scattering experiments so that the current PDG2010 [hi ] 
quotes \ V c d\ from neutrino charm production with an er- 
ror of about ~5%. With the new lattice calculations 
described in this article, the accuracy of \V c d\ from D 
semileptonic decays is approaching that from neutrino 
scattering and this provides an important consistency 
check. We find, 

\V cd \ =0.225(6)^.(10)^., (1) 

where the first error is from experiment [l2T ] and the 
second is the theory error from the Lattice QCD cal- 
culation presented here. Eq.([T]) is in excellent agree- 
ment with the PDG value based on neutrino scattering 
of \V cd \ = 0.230(11). 

In the rest of this article we describe how the result of 
eq.([T]) was obtained. We work with HISQ valence charm 
and light quarks on the MILC AsqTad Nf = 2 + 1 coarse 
(a ~ 0.12fm) and fine (a ~ 0.09fm) lattices. Table I lists 
the five MILC [lj| ensembles employed in this work and 
some simulation parameters. Compared to [1 01 ] we have 
doubled the statistics on ensembles CI, C2 and Fl. The 
valence charm and light bare quark masses are the same 
as in (lo| with the former tuned to the r) c mass and the 
latter chosen such that the ratio of mu g ht to the phys- 
ical strange quark mass is approximately the same for 
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TABLE I: The MILC N f = 2 + 1 ensembles used in the 
D — > 7r semileptonic project. The fourth column gives the 
valence HISQ light and charm quark masses in lattice units. 
N con f is the number of configurations and N tsr c the number 
of time sources used for each configuration. 



Set 


n/a 


mi(sea) /m s {sea) 


QiTTl valence 


N con f 


Ntsrc 


L 3 x N t 


CI 


2.647 


0.005/0.050 


0.0070 
0.6207 


1200 
1200 


2 
2 


24 3 x 64 


C2 


2.618 


0.010/0.050 


0.0123 
0.6300 


1200 
1200 


2 
2 


20 3 x 64 


C3 


2.644 


0.020/0.050 


0.0246 
0.6235 


600 
600 


2 
2 


20 3 x 64 


Fl 


3.699 


0.0062/0.031 


0.00674 
0.4130 


1200 
1200 


4 
4 


28 3 x 96 


F2 


3.712 


0.0124/0.031 


0.0135 
0.4120 


600 
600 


4 
4 


28 3 x 96 



valence and sea quarks. In the next section we summa- 
rize the formulas for hadronic matrix elements necessary 
to extract f+~**(0) and explain how they are related to 
three- and two-point correlators evaluated numerically 
on the lattice. These relations are the same as those 
described in reference [lOj so we will be brief. In sec- 
tion III we describe our data analysis and fitting proce- 
dures. We employ Bayesian fitting methods and carry 
out multi-exponential fits to several three-point and two- 
point correlators at the same time. Section IV discusses 
chiral and continuum extrapolations of lattice results to 
the physical limit. We apply the modified z-expansion 
method developed for D —> K semileptonic form factors 
in 10]. Section V presents our results for /+~ Hr (0) and 
| V cc i | and comparisons with other determinations of these 
quantities. Section VI gives a brief summary and we also 
include a 2 nd row unitarity test with all theory inputs 
coming from Lattice QCD. 



II. 



RELEVANT MATRIX ELEMENTS 



The most efficient way to calculate f+(q 2 ) at q 2 = 
is to focus on the scalar form factor fo(q 2 ) and use the 
kinematic identity /+(0) = /o(0). The scalar form factor 
can be determined from the matrix element of the scalar 
between the D meson and pion states. 



current S = * g 'F 



(m 0c - m i)(w\S\D) 



M 2 D 



Mi 



(2) 



The combination in the numerator of eq.([2]) does not get 
renormalized. The use of absolutely normalized currents 
is one of the reasons why we are able to significantly 
reduce errors in our D semileptonic scalar form factor 
calculations, both here and in [lOj . 



Our goal is to determine the hadronic matrix element 
(ir\S\D) ineq.([2]) via numerical simulations. The starting 
point is the three-point correlator, 



c^W^^EEE 



= ip x -(z-x) 



L 3 

x y z 

to) §{z,t) & D (y,t -T)}. 



(3) 



In eq.((3]) the interpolating operator $^ creates a D me- 
son at time slice to — T. At time t (to > t > to — T) the 
scalar current S converts the c quark inside the D into a 
light quark and also inserts momentum p w . The resulting 
pion then propagates to time slice to where it is annihi- 
lated by $ w . In addition to the three-point correlator one 
needs the pion and D meson two-point correlators, 



c 2 D pnt (t,t ) = ^eE^gm)*^*'*"))' (4) 



and 



C*> nt (t,t ;p«) = 

JJ E E e ^' iS ~ Vl h)). (5) 



Details on how the above three- and two-point correlators 
can be expressed in terms of single component staggered 
quark propagators are given in section IV of reference 
[10l | and will not be repeated here. There we also show 
how the sums J2 S in eqns.(|3]), (j4]) and (O can be carried 
out using U(l) random wall sources. 

The meson creation operators $|j and $J create not 
only the ground state D and pion we are interested in 
but also an arbitrary number of excited states with the 
same quantum numbers. Hence the t dependence of the 
two- and three-point correlators is complicated especially 
for staggered quarks. For two-point correlators it is given 
by, 

iV D -l 

C 2 D pnt (t) = £ bfie-^ + e-*?^-*)) 

3=0 
N'n-l 

+ E d^(-iy(e- E '^ + e- E '°^-% 

(6) 

and similarly for C 2pnt (t), except that there is no oppo- 
site parity terms for zero momentum. For three-point 
correlators one has 

N*—l Nn-l 

C 3 P nt (t,T)= J2 E Ajke-^e-^P-V 

N n -1 N' D -l 

+ E E B jk e- E Pe- E '* D V-*X-l)V-V 

j=0 k=0 
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+ E E C jk e-^ e - E ^(-iy 

j=0 k=0 
N'^-l N' D -1 

+ E E Djke-^e-^^i-lfi-l)^. 

j=0 k=0 

(7) 

We are interested in the ground state contributions with 
amplitudes, 



\^d\D)\ 2 
2M D a 3 



(8) 



(9) 
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FIG. 1: Ground state pion energy in lattice units for mo- 
mentum 2f (1, 1, 1) versus the number of exponentials in the 
fit. 



So the hadronic matrix element (7r| S"|Z?) that enters into 
the formula for fo'~ yiT (0) in @ is given by, 



(n\S\D) = 2y/M D E v 



D 



(11) 



We have accumulated simulation data for zero momen- 
tum D correlators and for pion correlators with momenta 
^(0,0,0), % (1,0,0), % (1,1,0) and f (1,1,1). The 
corresponding three-point correlators were calculated for 
T = 15, 16 on the coarse and for T = 19, 20, 23 on the 
fine ensembles. In the next section we describe how the 
combination on the right-hand-side of (jlip is obtained 
from three- and two-point correlators. 



III. 



FITS AND DATA ANALYSIS 



Extracting energies and amplitudes from numerical 
data on two- and three-point correlators is one of the 
more challenging but at the same time very important 
aspect of lattice calculations. For the past decade the 
HPQCD collaboration has been employing fitting meth- 
ods based on Bayesian statistics and involving multi- 
exponential fits [16| . For instance, in order to obtain the 
ground state energy and amplitude from a two-point cor- 
relator, we drop the first 1 ~ 4 time slices and then fit to 
the form of eq.© for several values of N = Nu (or N^) 
and N' . One continues to increase the number of expo- 
nentials until the fit results for Eq and feo including their 
errors and the chisquared per degree of freedom of the 
fit have stabilized. Fig.l shows an example of aE^ ver- 
sus N exp — N-n for momentum -J (1, 1, 1) with = N n . 
One sees that fits have stabilized after N exp = 3. It 
should be noted that as one increases the number of ex- 
ponentials and with it the number of fit parameters, the 
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FIG. 2: The form factor fo{q 2 ) for pion momentum ^(1, 1, 1) 
from simultaneous fits to two- and three-point correlators ver- 
sus N exp = Nd = N' D — N n — N^. The "sequential fitting" 
method was employed to go from one N exp value to the next. 



is inserted as initial conditions for the subsequent [N exp 
exponential fit. 



+ 11 



number of data points is growing as well. Each new fit 
parameter adds another prior term, i.e. additional data, 
to the fit function and the number of data points minus 
the number of fit parameters remains constant [l6l |. 
The combination „ is obtained from simultane- 

ous fits to C% nt {t), Cl pnt {t) and C 3pnt (t,T) for 2 (or 3) 
T values. In order to be able to include a large number 
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ID 



N exp and 



2 or N exp = 3. The 



of exponentials in these complicated fits we proceed as 

follows. We set N n = = N, 

start out with a small value, N exp 

fit results are then inserted as initial conditions for the 

subsequent N exp + 1 exponential fit. This procedure is 

repeated until one has completed multi-exponential fits 

2 shows results for 
2ttM 1 i\ ver _ 

One 



with N exp much larger than 2. Fig. 
foil 2 ) on ensemble C2 at pion momentum ^ (1,1,1 
sus N exp using this "sequential fitting" procedure, 
sees that, similar to in Fig.l, fit results have stabilized 
for N exp > 3 [H. 

In ongoing work we are investigating further methods 
to deal with complicated fits with large number of param- 
eters, in particular fits to collections of sums of exponen- 
tials For the calculations of this article, however, 
we have found that the "sequential fitting" method de- 
scribed above works well for all our data [llj]. We are 
even able to fit data on a given ensemble for all four 
pion momenta simultaneously and this allows us to ob- 
tain correlations between form factor results at different 
q 2 . These simultaneous fit results for foip-n) are given in 
Table II for several pion momenta p„ (the latter in units 
of 2*0. 



IV. CHIRAL AND CONTINUUM 
EXTRAPOLATION 

The next step is to extrapolate the data of Table II to 
the chiral/continuum limit. We do so using the "modified 
z-expansion fit" developed in 10]. The scalar form factor 
is parameterized as 



fo(q 2 ) 



(a D + a\D x z + a 2 D 2 z 2 



P(q 2 )$o 

x(l + b 1 (aE w ) 2 + b 2 (aE n ) 4 ), (12) 



with 



Di = 1 + c\xi + c 2 xilog(xi) + di(am c ) 2 + ei(am c ) 4 
'1 



+f i [l8M% + 8Mj c ), 



Ml 



xi = 



(13) 



(14) 



SM. 



5M 2 K = 



1 



{(M° ea ) 2 - (M™ /ence ) 2 ) , (15) 



((Mff Q ) 2 - [M v K alence ) 2 ) . (16) 



(47T/.) 2 

The kinematic variable z is defined as 20 2/ 



z(q 2 ,t Q ) = 



- q 2 - yt+ - h 

- q 2 + - t 1 



(17) 



TABLE II: Results for fo(p-w) for each ensemble 



Set 


fo (0,0,0) 


/o(l,0,0) 


fo (1,1,0) 


fo (1,1,1) 


CI 


1.1557(74) 


0.9155(93) 


0.8119(68) 


0.7384(209) 


C2 


1.1014(38) 


0.8700(59) 


0.7801(43) 


0.7315(87) 


C3 


1.0398(28) 


0.8787(36) 


0.7922(31) 


0.7326(66) 


Fl 


1.1053(29) 


0.8652(52) 


0.7586(53) 


0.7112(108) 


F2 


1.0443(28) 


0.8613(31) 


0.7645(53) 


0.7097(70) 



1.2 
1.1 
1 

,0.9 
0.8 
0.7 
0.6 



chiral/conti. extrapolation 
C3 ensemble 
C2 ensemble 
CI ensemble 

f atq 2 = 








0.2 



0.4 



0.6 
[GeV 2 



0.8 



FIG. 3: Chiral/continuum extrapolation of fo^ 71 versus E%. 
The data points are from the coarse ensembles (CI, C2 and 
C3). The three individual curves and the extrapolated band 
are from a fit to all live ensembles. 



with to a free parameter (which we set equal to 1.95GeV 2 
) and t± = (M D ± M^) 2 . We take $ from [H| and set 
P(q 2 ) = 1, where the latter relation reflects the absence 
of subthreshold poles in the scalar channel. 

We show results of fits to the form of eq. (fT2"|) in Figs. 3 
& 4. We plot separately the coarse and fine data points in 
order to be able to better distinguish individual curves. 
However, the fit was done simultaneously to all the data 
in Table II, coarse and fine. The y^/dof = 0.85 for 
this fit. In Fig. 5 we show both fo^ K {q 2 ) (results from 
[101 ) and fo^^iq 2 ) versus q 2 in the physical region 
< q 2 < q 2 nax — [M D — M 7r ( A -)) 2 . Note that we consider 
the correlations between the momenta in the fits which 
we did not consider in |10| . However, we find that includ- 
ing or excluding correlations in our chiral/continuum ex- 
trapolations has minimal effect on /J 3 ^ x (0) at the phys- 
ical point, namely a ~ 0.04<j shift in the central value 
and a ~ 0.02c change in the error. 

The motivation for the "modified z-expansion fit" is ex- 
plained in more detail in [lOj. Form factors at q 2 = are 
needed to extract the CKM matrix elements \ V c d \ or |V CS |. 
The pion energy approaches lGeV in this kinematic re- 
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FIG. 4: Same as for Fig. 3 showing, however, the fine data 
points. 



TABLE III: Error Budget for /+~ Mr (0) 



Type 


Error (%) 


Statistical 


Q 1 

6.1 


Scale (ri and T\j<x) 


U. / 


Expansion coeff. d% 


U.o 


FT 


u.o 


i^rgiiL quarK. Liepeiioence 


1 Q 


Sea quark dependence 


0.6 


am c corrections 


2.0 


aE w corrections 


1.0 


Finite volume 


0.04 


Charm mass tuning 


0.05 


Total 


4.4% 



physical limit consistent with the z-expansion result and 
with comparable errors, however with worse % 2 / dof . 




0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 
q 2 [GeV 2 ] 



FIG. 5: fo-+ K {q 2 ) and fa^(q 2 ) versus q 2 in the physical 
limit. 



gion and so chiral perturbation theory might cease to be 
valid. The z-expansion, on the other hand, is applica- 
ble throughout the physical kinematic region, and our 
new "modified z-expansion" allows for the expansion co- 
efficients to be mass and lattice spacing dependent. We 
have checked that fits to eq. (fT2|) are stable with respect to 
adding further terms such clS ■ (am c ) 6 , (oE'tt) 6 or keep- 
ing just the cj and c\ terms in eq. (fl"3)) (the ZVs). Such 
changes in the fit ansatz led to minimal changes in both 
the central value and the error for /n(0) in the physical 
limit. We have also verified that traditional ChPT ex- 
trapolations (see Appendix C & D of [l(| and references 
therein for relevant ChPT formulas) lead to /+(0) in the 



V. RESULTS FOR f+~* n (Q) AND \V cd \ IN THE 
PHYSICAL LIMIT 

Our final result for the D —> ir form factor at q 2 = 
averaged over D° ir~ and D + —> ir°, is, 

/^(0) = 0.666(20) stat .(21) sys ., (18) 

where the first error is statistical and the second sys- 
tematic. Fig. 6 plots our new result together with other 
theory calculations [3, US H3] and experimental deter- 
minations [H, HI] (the latter use CKM unitarity values 
for \V cd \ to extract f+-**(0)). 

The total error in our /r )_> " 7r (0) is 4.4%, signifying a 
better than factor of 2 improvement over previous lattice 
determinations. The full error budget is given in Table 
III. The largest error is statistical followed by am c and 
light quark mass dependence errors. All but the last 
two entries in Table III were obtained using the methods 
described in reference [26[ and Appendix B of For 
instance, the "light quark dependence" errors come from 
the c\ and c\ terms in the fit ansatz eq.(13), the "am c 
corrections" from the di and ej terms etc. 

Finite volume errors were estimated by calculating a 
pion tadpole integral both at finite and at infinite volume. 
The charm mass tuning error is taken to be the same as 
for D — > K, Iv for which calculations at two values of am c 
were carried out explicitely to estimate this error. Effects 
from electromagnetism/isospin breaking and charm sea 
are expected to give negligible contribution to the error 
budget compared to our other errors (see (Io||). 

Finally one can combine our result for /P~*' 7r (0) with 
CLEO-c's measurement of \V cd \ x /+ ^(O) to extract 
a precision value for \V c d\ from D scmilcptonic decays. 
This leads to the result quoted already in eq.([T]), which 
is shown in Fig. 7 together with values from PDG2010. 
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HPQCD (2011) 
f ^(O) = 0.666 (20) (21) 

Fermilab/MILC (2005) 
Sum Rules (2009) 
ETMC (2011, preliminary) 

Experiment + CKM Unitarity 
CLEO-c (2009) 

Belle (2006) 
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4 % error 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
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FIG. 6: The D -> vr form factor /£ ^(OWrom this work and 
comparisons with other determinations [l^ . Il3l . I23r - t25| . 



HPQCD 2011 
IV J = 0.225 (6) exp (10), aL 



PDG 2010: Semileptonic decay h 
(Fermilab/MILC 2005) 

PDG 2010: Neutrino exp. 

PDG 2010: Unitarity 



0.1 



0.15 



0.2 



cd 



0.25 



FIG. 7: Comparison of \V c d\ from this work with values in 
PDG2010 0. 



VI. SUMMARY 

In this article we have presented a new calculation of 
the D — > tt,Il> semileptonic form factor f+^^iq 2 ) at 



q 2 = 0, with errors a factor of two better than in the 
past. This combined with recent precision measurement 
of the branching fraction for this process by CLEO-c has 
allowed for an accurate determination of the CKM ma- 
trix element \V c d\. Direct determination of \V c d\ from D 
semileptonic decays is becoming competitive with that 
from neutrino scattering. The fact that these two very 
different processes lead to the same \V cc i\ is a nontrivial 
consistency check of the Standard Model. 

Finally using our values for \V c d\ and |V CS | plus 
the most recent |T4j,| e:cc j. = 0. 0397(10) from the Fermilab 
Lattice & MILC collaborations [231, the most up-to-date 
test of second row unitary from Lattice QCD becomes, 



W, 



cd 



Wr 



\ 2 + \v cb \ 2 



0.976(50). 



(19) 



This improves on the PDG2010 value 1.101(74) [Tj]. 

In the future we will be reducing the largest errors 
in Table III by increasing statistics and simulating on 
finer lattices [28| . Calculations of the full q 2 dependence 
of f+^ K (q 2 ) and f+^(q 2 ) are also already underway 
28). Furthermore we are working on updating HPQCD's 
result for the D meson decay constant fn || and on 
carrying out an independent extraction of \V c d\ from D 
leptonic decays (29j . 
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Appendix A: Priors and Prior Widths for Two- and 
Three-Point Correlators 



In this appendix we give sample priors and prior widths 
used in the fits of section III (the reader is referred to 
reference [l6| for definitions of these terms). We turn 
to the Particle Data Group listings for guidance in pick- 
ing priors for energies and energy splittings. All energies 
in Table IV are given in lattices units and are appro- 
priate for ensemble C2. Numbers for other ensembles 
can be obtained by rescaling with corresponding lattice 
spacings. Prior widths for amplitudes are fixed based on 
exploratory initial fits. 
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